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In this section is presented a new integer number algorithm for linear equation. 
This algorithm is more “rapid” than W. Sierpinski’s presented in [1] in the sense that it 
reaches the general solution after a smaller number of iterations. Its correctness will be 
thoroughly demonstrated. 

INTEGER NUMBER ALGORITHM TO SOLVE LINEAR EQUATIONS 



Let’s us consider the equation (1); (the case a^b e©, i = \,n is reduced to the 
case (1) by reducing to the same denominator and eliminating the denominators). Let 
d = (a ] ,...,a n ) .If d\b then the equation does not have integer solutions, while if d\b 
the equation has integer solutions (according to a well-known theorem from the number 
theory). 

If the equation has solutions and d & we divide the equation by d . Then d = 1 
(we do not make any restriction if we consider the maximal co-divisor positive). 

Also, 

(a) If all a ] the equation is trivial; it has the general integer solution 



x,. = k t eZ, i = 1 ,n, when b = 0 (the only case when the general solution is 
n -times undetermined) and does not have solution when/? ^ 0 . 

(b) If 3i, 1 < i < n such that a, = ±1 then the general integer solution is: 

r \ 



x , - -a, 



ajkj 



\j*i 



and x s = k, e LI , s e 1, 



n \ i 



The proof of this assertion was given in [4]. All these cases are trivial, therefore 
we will leave them aside. The following algorithm can be written: 



Input 

A linear equation: 

n 

(2) ^ ] a i x j = b, a n b eZ, ^ ±1, i = 1 ,n , 

Z = 1 

with not all a t = 0 and (a ] ,...,a n ) = 1 . 
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Output 

The general solution of the equation 

Method 

1. h := 1, p := 1 

2. Calculate min Irl, r = a,, (mod a,), Irlca. and determine r and the pair 

1 <i,j<n J ' J ' 

(. i,j ) for which this minimum can be obtained (when there are more possibilities we have 
to choose one of them). 

3. If \r\ ^ go to step 4. 

If Irl = 1 , then 



-V- := r 



//, - X + b 



v 

f 



s = 1 
s<£ i,j 



Cl- — V A — U- 

a i t h+ X a s X s+ b 



\Xj:=r\ 



a j S=1 

ij 



r -a. 
a. 



j 



(A) Substitute the values thus determined of these unknowns in all the 
statements ( p ), p = 1,2,... (if possible). 

(B) From the last relation ( p ) obtained in the algorithm substitute in all 
relations: (p - l),(p - 2),...,(1) 

(C) Every statement, starting in order from ( p - 1) should be applied the 
same procedure as in (B): then {p - 2),. ..,(3) respectively. 

(D) Write the values of the unknowns jq, i = l,n, from the initial 
equation (writing the corresponding integer number parameters from 
the right term of these unknowns with k x ,...,k n _ x ), STOP. 

a — r 

4. Write the statement ( p ) : x. = t h — x, 

a j 

5. Assign Xj'.= t h h:=h + 1 

a i := r p := p + \ 

The other coefficients and variables remain unchanged go back to step 2. 

The Correctness of the Algorithm 

Let us consider linear equation (2). Under these conditions, the following 
properties exist: 
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Lemma 1. The set M = |r| , r = a j (mod a j ), 0< |r| < | | has a minimum. 

Proof: 

Obviously McN’ and M is finite because the equation has a finite number of 
coefficients: n , and considering all the possible combinations of these, by twos, there is 
the maximum AR 2 (arranged with repetition) = n 2 elements. 

Let us show, by reduction ad absurdum, that M ^ 0 . 

M ^ 0 o a, = 0(moda y ) Vz, j = 1, n. Hence a . = 0(mod a t ) Vz, j = 1 ,n. Or this 
is possible only when \a\ = V/, j = \,n , which is equivalent to 

(a l ,..,a n ) = cij, Vi el,n . But (a l ,..,a n ) = l are a restriction from the assumption. It 
follows that \a\ = \,n, Vi el ,n a fact which contradicts the other restrictions of the 
assumption. 

M ^ 0 and finite, it follows that M has a minimum. 

Lemma 2. If \r\ = min M , then Irl < la, I, Vi el, a . 

1 <i,j<n 

Proof: 

We assume conversely, that 3i 0 , 1 <i 0 <n such that |r| > | a io | . 

Then |r| > min a. = a , Ut 1 , 1 < / n < n . Let a„ , 1 <p n <n, such that \a„ > a, and 

1<;'<« I 1 1 | A | Po r u I Pa I I | 

a ■ is not divided by a 0 . 

Po J J 

There is a coefficient in the equation, |a y | which is the minimum and the coefficients 

are not equal among themselves (conversely, it would mean that (aj,..,a„) = a x = ±1 
which is against the hypothesis and, again, of the coefficients whose absolute value is 
greater that |a i)o | not all can be divided by a ; (conversely, it would similarly result in 
(a l ,..,aJ = a /o * ±1 . 

We write [a po / a /o j = q t) e Z (integer portion), and r = a po - cy 0 a ; e Z . We have 

a /)( = r () (mod a /( ) and 0 < |r 0 | < | a k | < |a | < |r| . Thus, we have found an r 0 which |r 0 | < |r| 
contradicts the definition of minimum given to |r| . 

Thus|r| <\a\, Vi el ,n . 

Lemma 3. If |r| = minM = 1 for the pair of indices (/,/) , then: 
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x, = r 



- E fl A +h 



Xj = r 



5^=1 

i,j 



\ 



a.-r ^ r-a. 

a,t h + E fl A + /7 

a. , =1 a. 

$£ 1,7 



x = k eZ, se l,...,n \ i, j 



is the general integer solution of equation (2). 

Proof: 

Let x =x° e=l,n, be a particular integer solution of equation (2). Then 



y 

3k s = x s ° eZ, s e {l,...,«}\ {/,y } and t h = x“ + — x,° eZ (because a,. - r = Ma y . ) such 

that: 



a y 



x, = r-a. 



Xj = r — a j 



x j+- 



a. - r 



CI J j 



x j+- 



a - r 



y 



- E a A° +/:; = T 

n 

-z 



5=1 

i,j 



+ 



a. - r 



a j 5=1 

•se i,j 



0 , r-a, J 0 

a s x s + = x, 

a. 



and x ( = fc = x° , s e \ i, j . 



Lemma 4. Let \r\ ^ and (/, /) be the pair of indices for which this minimum can 
be obtained. Again, let’s consider the system of linear equations: 



(3) 



a i t h +rx i+ E a s X s = b 

5=1 
5£ i,j 

a. -r 

t h= x j+ x i 

a j 



Then x e -x ) e , e = l, n is a particular integer solution for (2) if and only if x e = x° e , 
o 0 a. — r 

e e 1 \ j and t h = t h = x . + — x. is the particular integer solution of (3). 

“j 

Proof: 

x e =x° e , e = l,n is a particular solution for (2) if and only if 
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'Ea e x° e =bd £ a s x° s + 

e=l s = 1 

^ i, 7 

n 

<=>«/" + ^,° + Z a v X “ =/? 



x ,+- 



a. - r 



a ; 



+ rx° = b o 



and 



ti= x ° j+ ^: x ° eZ 



<=> x e = x , 



«=i 
*£ i,j 



a 



j 



ee 1 ,...,« \ j and t h = is a particular integer solution for (3). 



Lemma 5. The previous algorithm is finite. 

Proof: 

When \r\ = 1 the algorithm stops at step 3. We will discuss the case when 
\r\ ^ 1 . According to the definition of r, |r| eN* . We will show that the row of 
r — s successively obtained by following the algorithm several times is 
decreasing with cycle, and each cycle is not equal to the previous, by 1. Let r t be 

the first obtained by following the algorithm one time, q | ^ 1 then go to step 4, 

and then step 5. According to lemma 2, r, | < \a\, \/i = 1 ,n . 

Now we shall follow the algorithm a second time, but this time for an 
equation in which q (according to step 5) is substituted by a, . Again, according 

to lemma 2, the new |r| written |r 2 | will have the propriety: \r 2 \ < q . We will get 
to \r\ = 1 because \r\ > 1 and \r\ < oo, and if \r\ ^ 1 , following the algorithm once 
again we get |r|<|q| and so on. Hence, the algorithm has a finite number of 
repetitions. 

Theorem of Correctness. The previous algorithm calculates the general 
solution of the linear equation correctly (2). 

Proof: 

According to lemma 5 the algorithm is finite. From lemma 1 it follows 
that the set M has a minimum, hence step 2 of the algorithm has meaning. When 
\r\ = 1 it was shown in lemma 3 that step 3 of the algorithm calculates the general 
integer solution of the respective equation correctly the equation that appears at 
step 3). In lemma 4 it is shown that if \r\ ^ 1 the substitutions steps 4 and 5 
introduced in the initial equation, the general integer solution remains unchanged. 
That is, we pass from the initial equation to a linear system having the same 
general solution as the initial equation. The variable h is a counter of the newly 
introduced variables, which are used to successively decompose the system in 
systems of two linear equations. The variable p is a counter of the substitutions of 
variables (the relations, at a given moment between certain variables). 

When the initial equation was decomposed to \r\ = 1 , we had to proceed in 
the reverse way, i.e. to compose its general integer solution. This reverse way is 
directed by the sub-steps 3(A), 3(B) and 3(C). The sub-step 3(D) has only an 
aesthetic role, i.e., to have the general solution under the form: x t = j]{k [ ,...,k n _ x ) , 
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i = \,n, f] being linear functions with integer number of coefficients. This “if 
possible” shows that substitutions are not always possible. But when they are we 
must make all possible substitutions. 

Note 1. The previous algorithm can be easily introduced into a computer 
program. 

Note 2. The previous algorithm is more “rapid” than that of W. 
Sierpinski’s [1], i.e., the general integer solution is reached after a smaller number 
of iterations (or, at least, the same) for any linear equation (2). 

In the first place, both methods aim at obtaining the coefficient +1 for at least 
one unknown variable. While Sierpinski started only by chance, decomposing the 
greatest coefficient in the module (writing it under the form of a sum between a 
multiple of the following smaller coefficient (in the module) and the rest), in our 
algorithm this decomposition is not accidental but always seeks the smallest \r\ 
and also choose the coefficients a, and a . for which this minimum is achieved. 

That is, we test from the beginning the shortest way to the general integer 
solution. Sierpinski does not attempt to find the shortest way; he knows that his 
method will take him to the general integer solution of the equation and is not 
interested in how long it will take. However, when an algorithm is introduced into 
a computer program it is preferable that the process time should be as short as 
possible. 

Example 1. 

Let us solve in Z 3 the equation \lx — ly + lOz = -12 . 

We apply the former algorithm. 

1. h = 1 ,p = 1 

2. r = 3,i = 3,j = 2 

3. |3| ^ 1 go on to step 4. 

„ /1X 10-3 

4- (1) y = t x —■z = t l + z 

5. Assign 

y := t x h := 2 

a 3 := 3 p:= 2 

with the other coefficients and variables remaining unchanged, go back to 

step 2. 

2. r = — 1, i' = l, j =3 

3. |-1| = 1 

x = -l(-3 1 2 - 12) = 3 1 2 - 7fj - 12 

z = -l( \lt 2 + (-70 • 17 ~ ( ~ 1} + =^(-12)] = 11 1 2 + Alt, - 72 
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(A) We substitute the values of x and z thus determined into the only 
statement ( p ) we have: 

(1) y = t x +z = — 17^+43^-72 

(B) The substitution is not possible. 

(C) The substitution is not possible. 

(D) The general integer solution of the equation is: 

x = 3k x - lk 2 + 12 
< j = -17^+43^ -72 
z--Ylk x +42^ -72; k x ,k 2 eZ 



REFERENCES: 

[1] Sierpinski, W, - Ce stim si ce nu stim despre numerele prime? - Editura 
Stiintifica, Bucharest, 1966. 

[2] Creanga, I., Cazacu, C., Mihut, P., Opait, Gh., Corina Reisher - 
Introducere in teoria numerelor, Ed. Did. si Ped., Bucharest, 1965. 

[3] Popovici, C. P. - Aritmetica si teoria numerelor, Ed. Did. si Ped., 
Bucharest, 1963. 

[4] Smarandache, Florentin - Un algoritm de rezolvare in numere intregi a 
ecuatiilor liniare. 



7 




